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ABSTRACT 


The  results  of  a  preliminary  study  on  thermal  modeling  and  analysis 
of  microcircuits  are  presented.  Studies  covering  typical  components 
mounting  methods,  bonding  agent  materials  component  spacing,  and 
cooling  concepts  are  described.  Temperature  profiles  of  substrates, 
component  temperatures  and  thermal  resistances  are  presented.  It 
is  shown  that  the  use  of  moly  tabs,  sputtered  beryllia  and  thermo¬ 
electric  coolers  can  significantly  reduce  component  temperature 
levels.  Additionally,  it  is  illustrated  that  thermal  modeling 
can  aid  the  circuit  designer  in  the  layout  and  preliminary  design 
phases. , 
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1.0  INTRODUCTION 

Microcircuits  consist  of  integrated  circuit  devices,  semi-conductor 
components,  and  other  non- semiconductor  elements  mounted  in  a 
patterned  array.  The  various  devices  are  electrically  inter¬ 
connected  by  thick  film  conductors,  metalized  connecting  strips  or 
wire  bonds.  The  array  is  attached  to  a  suitable  substrate  such 
as  alumina,  beryllia,  or  porcelain.  Typically  a  substrate  is  then 
bonded  to  a  base  material  and  placed  in  a  hermetically  sealed 
container,  as  shown  in  Figure  (1). 

Primary  sources  of  heat  generation  are  considered  to  occur  at 
emitter  junction  regions  of  the  semiconductor  devices.  Secondary 
sources  for  thermal  dissipation  occur  in  the  discrete  elements 
such  as  thick  and  thin  film  resistors.  Since  both  large  temperature] 
excursions  and  thermal  cycling  can  significantly  decrease  the  life 
of  the  circuit,  consideration  of  thermal  effects  must  be  included 
during  the  design  phase. 

The  importance  of  thermal  cycles  and  increased  circuit  temperatures 
may  be  illustrated  by  considering  life  testing  techniques  for  semi¬ 
conductors  (1)*.  Accelerated  life  testing  can  be  achieved  by 
thermal  cycling.  That  is,  through  thermal  cycling  the  life  time 
span  of  normal  equipment  can  be  compressed  by  a  factor  of  200  to 
400.  Similarly,  elevated  temperatures  can  result  in  increased 
aging.  Silicon  devices  operated  at  300°C  will  age  10,000  times 
faster  than  tt-*>  same  device  at  125°C.  The  effect  of  a  200°C 
environment  compared  to  one  at  125°C  results  in  increased  aging  at 
a  rate  of  170  times  faster.  Thus,  the  minimization  of  both 
temperature  gradients  and  thermal  cycling  should  be  a  constraint 
of  the  circuit  designer. 

The  circuits  considered  herein  are  usually  mounted  on  a  substrate 
by  either  a  metal  filled  epoxy  or  by  a  eutectic  bond.  Transfer 
of  the  thermal  energy  dissipated  by  the  components  is  primarily  by 
conduction  through  the  bond  material  into  the  substrate.  The 
thermal  energy  then  diffuses  through  the  remaining  bond  material 
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into  the  base  where  the  heat  can  be  removed  from  the  package.  Tht 
removal  of  the  thermal  energy  from  the  package  base  can  be 
acheived  by  a  variety  of  standard  techniques  constrained  only  by 
the  overall  packaging  method.  That  is,  a  conductive  path  through 
a  circuit  card  to  the  mounting  structure  can  be  employed  or  forced 
convection  exchange  with  the  ambient  environment  may  be  utilized. 

Circuit  components  are  extremely  small  in  physical  dimensions  but 
their  power  dissipation  can  be  relatively  large  in  terms  of  heat 
flux.  It  is  not  uncommon  for  a  "power  chip"  to  have  dimensions 

of  20  mils  by  15  mils  while  dissipating  nearly  one  watt  of  energy. 

6  2 

This  gives  a  heat  flux  of  1.64  X  10  BTU/hr  ft  over  an  area  of 
0.0003  square  inches.  The  combination  of  the  large  localized 
heat  fluxes,  state  of  the  art  mounting  techniques,  and  proposed 
large  scale  integrated  circuit  (high  density  packaging)  are  not 
compatible  with  current  cooling  techniques.  Thus  some  technique 
is  needed  by  the  designer  to  determine  the  best  spatial  locations 
of  the  components  in  the  array  and  for  investigating  new  mounting 
and  cooling  methods.  Thermal  modeling  offers  such  as  design 
tool  permitting  parametric  design  changes,  temperature  predictions 
and  comparative  temperature  control  studies  to  be  accomplished 
prior  to  circuit  fabrication  (2),  (3),  (4),  and  (5). 
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2.0  ANALYSIS 

A  physical  description  of  a  typical  hybrid  microcircuit  and 
packaging  container  is  illustrated  by  the  sketch  given  in 
Figure  1.  Each  circuit  contains  multiple  semiconductor  devices 
and  resistor  elements  mounted  on  a  common  base  which  may  be  cooled 
by  a  variety  of  techniques.  The  thermal  energy  dissipated  by  the 
multiple  components  and  the  presence  of  various  material  and 
interfaces  presents  a  heat  transfer  problem  which  can  be  solved 
realistically  only  by  computer  methods,  (6)  and  (7). 

A  generalized  lumped  parameter  nodal  thermal  model  was  formulated 
to  represent  the  various  materials,  interfaces  and  multiple  heat 
dissipating  components  of  a  hybrid  microcircuit,  (6).  The  physical 
system  was  sectioned  into  an  array  of  discrete  nodal  volumes,  and 
the  mass  of  each  volume  was  then  "lumped"  at  a  point  within  the 
volume  which  it  represented.  The  paths  for  heat  transfer  from 
one  node  to  another  were  represented  by  conductors  joining  the 
appropriate  nodes.  The  thermal  model  contained  2000  nodes  to 
represent  the  circuit  and  substrate  layers.  The  nodes  were 
divided  into  8  layers  of  material  (volume  containing)  nodes  and  2 
layers  of  surface  (boundary)  nodes.  All  the  layers  were  arranged 
in  a  20  X  10  grid  pattern  with  the  size  of  the  divisions  in  the 
grid,  the  thermal  conductivity,  and  the  thickness  of  the  layers 
of  material  nodes  being  determined  by  the  program  user.  The 
lowest  layer  of  nodes,  which  represented  the  bottom  surface  of 
the  base,  were  boundary  nodes.  Above  the  layer  of  boundary  nodes 
were  four  layers  of  material  nodes,  thus  allowing  the  use  of 
four  layers  of  material  in  the  substrate.  Another  layer  of 
boundary  nodes  was  provided  on  the  top  surface  of  the  substrate. 
While  the  geometry  of  the  lower  six  layers,  or  substrate,  could 
only  be  that  of  a  square  or  rectangle,  it  was  possible,  through 
the  use  of  special  control  cards,  to  represent  various  shapes  and 
sizes  of  components  in  the  upper  four  layers  of  nodes.  These 
components  consisted  of  epoxy  layers,  powered  electrical  devices, 
resistors,  rooly  tab  mounts,  or  connecting  strips.  In  Figure  2 
an  example  circuit  is  shown  in  which  only  two  of  the  upper  four 
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3.0  DISCUSSION  OF  RESULTS 

A  parameter  which  is  often  specified  for  circuit  design  purposes 
is  the  thermal  resistance  of  a  mounted  powered  component.  The 
overall  thermal  resistance  is  defined  by  dividing  the  temperature 
difference  between  the  device  and  the  case  by  the  thermal  power 
dissipated  in  the  device.  In  the  case  of  a  one  watt  thermal 
dissipation  the  value  of  temperature  is  equivalent  to  the  thermal 
resistance.  Table  1  lists  typical  material  thermal  conductivities. 

3.1  Alumina  Substrates 

Using  the  analysis  methods  described  in  the  previous  section 
component  temperatures,  substrate  temperature  profiles  and  overall 
thermal  resistances  were  calculated  and  studied  for  a  variety  of 
cases.  The  first  case  investigated  concerned  typical  bonding 
methods  for  single  chips  to  alumina  substrates  with  an  epoxy 
bonding  agent,  as  shown  in  Figure  2.  This  particular  chip, 

20  X  30  X  5  mils  in  size,  was  considered  to  have  100  milliwatts 
of  thermal  energy  to  be  dissipated  through  the  laminate  sublayers. 

As  may  be  noted  in  Figure  3,  the  overall  thermal  resistance  between 
the  chip  and  kovar  base  is  greatly  influenced  by  the  size  of  the 
epoxy  bond.  Initially,  as  the  epoxy  layer  extends  beyond  the  edge 
of  the  chip,  the  thermal  resistance  decreases  until  the  fin 
effectiveness  of  the  epoxy  becomes  insignificant.  This  indicates  tha 
epoxy  "squeeze  out"  is  beneficial  until  the  area  of  the  bond  layer 
is  approximately  4  times  greater  than  the  chip  size  for  the  con¬ 
figuration  studied.  When  the  epoxy  area  equals  the  chip  area 
the  thermal  resistance  was  120°C  /  watt  while  at  an  area  ratio  of 
4  the  thermal  resistance  was  found  to  be  85°C/watt.  In  each 
specific  case  the  change  in  thermal  resistance  depends  upon  the 
ability  of  the  bonding  agent  to  diffuse  the  thermal  energy  over  a 
greater  area  prior  to  conduction  into  the  substrate.  It  is  this 
sensitivity  of  the  thermal  resistance  to  squeeze  out  which  some¬ 
times  causes  quality  and  temperature  reliability  problems.  That  is, 
circuits  which  are  thought  to  be  identical  will  have  components 
operating  at  different  temperatures. 
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MATERIAL 

THERMAL  CONDUCTIVITY* 
WATT/MIL-*C** 

SILICON 

0,005 

GOLD-SILICON  EUTECTIC 

0,0072 

ABELFILM517  CA1  EPOXY) 

& 00008 

OHMEX  SILVER  EPOXY 

0.0011 

GOLD  FILLED  EPOXY  (All  EPO,  Y) 

0,0001 

ALUMINA 

(100047 

ECCOBOND  99  EPOXY 

0.000043 

KOVAR 

aooo44 

MOLYBDENUM 

Q.0Q348 

•  VALUES  GIVEN  ARE  FOR  100*C 
~  I  WATT/MIL  *C  •  435  x  lfl4  BTU/HRFT  *F 


TABLE  1:  MATERIAL  THERMAL  CONDUCTIVITIES 
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layers  of  nodes  were  utilized.  The  upper  of  the  two  layers  con¬ 
tained  a  single  silicon  chip,  while  the  lower  layer  contained  an 
Ohmex-Ag  epoxy  pad. 

The  two  layers  of  surface  nodes  mentioned  in  the  previous  paragraph 
were  used  to  simulate  various  boundary  conditions.  That  is,  by 
placing  nodes  on  the  outer  surface  convection  and  radiation  boundary 
conditions  could  be  included.  Additionally,  the  model  allowed  any 
number  of  nodes  to  be  designated  as  constant  temperature  nodes, 
thus  fixing  their  temperature. 

Several  assumptions,  other  than  those  usually  required  for  nodal 
thermal  modeling,  were  made  during  the  studies.  It  was  assumed 
that  the  thermal  energy  of  the  various  components  was  dissipated 
uniformly.  Also,  thermal  energy  dissipated  by  or  conducted  through 
the  connecting  strips  was  of  secondary  importance.  Additionally, 
natural  convection  between  the  case  and  the  upper  substrate  surface 
was  negligible.  The  first  assumption  could  not  be  validated, 
however,  temperature  errors  introduced  by  this  assumption  would 
be  small  due  to  the  relatively  high  thermal  conductivity  of  the 
silicon  chips  in  comparison  to  the  substrate  ceramic.  A  comparison 
of  the  thermal  conductance  of  metalized  strips  to  that  of  the 
adjacent  substrate  was  made  and  the  metalized  strips  were  several 
orders  of  magnitude  less  than  through  the  alumina.  Natural 
convection  of  dry  nitrogen  in  the  small  cavity  (1  square  inch  by 
0.1  inch!  was  also  found  to  be  negligible.  For  this  case  the 
Grashof  number  was  31.1  and  is  too  low  to  be  significant  for  internal 
convection  effects. 

Preliminary  studies  investigating  the  utility  of  thermoelectric 
cooling  concepts  required  modification  in  the  thermal  model.  That 
is,  a  subroutine  which  incorporated  thermoelectric  fundamentals 
was  written  and  applied  in  terms  of  sources  and  sinks  at  node 
levels  which  normally  represented  the  bonding  agent  between  the 
substrate  and  Kovar  base,  Figure  (2).  Appendix  A  contains  a  descrip¬ 
tion  of  thermoelectric  fundamentals  and  a  listing  of  the  sub¬ 
routine  . 
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Typical  temperature  profiles  for  a  chip  bonded  to  a  laminate  are 
illustrated  in  Figure  4  and  5.  Figure  4  shows  the  importance 
of  minimizing  the  thermal  resistance  of  the  bonding  agent  between 
the  alumina  and  kovar  layers.  It  may  be  noted  that  the  aluminum 
filled  epoxy  in  this  case  tends  to  act  as  an  insulating  layer. 

Figure  5  shows  how  the  thermal  energy  diffuses  over  the  top 
surface  of  the  alumina.  In  both  figures  isotherms  are  shown  to 
depict  the  effective  thermal  paths. 

Figure  6  shows  the  effect  on  base  to  chip  temperature  difference 
for  two  identical  chips  mounted  on  the  same  substrate.  For  the 
case  under  consideration  the  chip  temperature  or  thermal  resistance 
decreased  as  the  chips  were  separated  by  increasingly  greater 
distances.  When  the  separation  was  approximately  equal  to  the 
bonding  pad  dimension  the  temperature  decreased  to  nearly  the 
value  of  a  singularly  mounted  component.  This  result  is  greatly 
influenced  by  the  bonding  pad  size.  Larger  bonding  pads  produce 
a  smaller  increase  in  temperature  when  mounted  close  together 
because  of  their  ability  to  diffuse  the  thermal  energy.  In  the 
case  illustrated  here,  a  pad  of  Ohmex-Ag  bonding  agent  50  mils 
X  50  mils  would  cause  only  a  3°C  increase  in  temperature  for  the 
same  spacing  criteria  in  comparison  to  the  7°C  found  for  the  30  x  30 
mil  bonding  pads. 

3.2  Alumina  Substrates  with  Heat  Conducting  Fluids 

The  above  results  indicated  that  a  potentially  attractive  method 
for  heat  removal  might  be  to  replace  the  inert  gas  which  covers 
the  circuit  face  with  a  heat  conducting  fluid.  Results  for  this 
study  are  given  in  Figure  (7).  As  may  be  noted  for  typical  ' 
silicone  oil  conductivities  of  0.5  BTU/hr  ft  °F  the  difference 
between  the  chip  temperature  and  the  base  was  reduced  by  approximate¬ 
ly  10%.  This  was  not  considered  to  be  a  significant  reduction  to 
warrant  the  added  complexity  introduced  by  such  a  solution. 

3.3  Moly  Tabs 

A  further  study  was  conducted  to  determine  the  temperature  control 
aspects  of  mounting  semiconductor  devices  on  moly  tabs.*  This 

|  *Moly  tabs  denote  Molybdenum  mounting  pads. 
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•  CHIP  TYPE  -  SILICON 

SIZE  -  20  X  20  X  5  MILS 
POWER  -  1  WATT 

•  EPOXY  PAD  TYPE  -  OHMEX-Ag 

SIZE  30  X  30  X  1  MILS 

•  SUBSTRATE  TYPE  -  ALUMINA-ECC080ND  99- ICO  VAR 

SIZE  -  200  X  300  X  50  MILS 


...  A  -J - 1 - 1 _ I 

20  30  40  50  60 

S-EPOXY  PAD  SEPARATION  —  MILS 


FIGURE  6:  EFFECT  ON  CHIP  TO  BASE  TEMPERATURE  DIFFERENCE  DUE  TO  SPACING  BETWEEN  CHIPS 
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mounting  technique  which  is  illustrated  in  Figure  8  provides  an 
effective  reduction  in  the  localized  heat  flux  similar  to  that 
achieved  by  the  epoxy  squeeze  out/  only  in  this  case  a  higher 
thermal  conductivity  material  which  can  be  pre-sized  is  utilized. 

The  results  for. the  thermal  resistance  of  the  various  size  chips 
are  given  in  Figure  f.  As  may  be  noted  the  moly  tab  tends  to 
diffuse  the  thermal  energy  in  a  manner  similar  to  that  of  a  fin. 
Initially  there  is  a  large  effect,  but  this  decreases  with  an 
increase  in  the  ratio  of  the  moly  tab  area  to  chip  area.  The 
computer  results  were  compared  with  a  measured  value  which  was 
experimentally  obtained  with  a  Bendix  Thermal  Analyzer.  The 
Bendix  Thermal  Analyzer  indicated  a  thermal  resistance  of 
approximately  56°C/watt  while  the  numerical  results  gave  a  value 
of  51°C/watt.  The  difference  between  the  two  values  may  be  due 
to  errors  in  physical  measurements,  inaccuracies  in  thermal 
conductivities,  or  due  to  an  unknown  thermal  resistance  at  a 
bond  interface. 

A  multiple  component  design  study  of  a  specific  circuit  was 
undertaken  to  investigate  various  component  mounting  methods.  The 
circuit  which  was  analyzed  is  shown  in  Figure  10,  and  a  description 
of  the  variations  considered  is  listed  in  Table  2.  The  component 
to  case  temperature  differences  for  the  four  cases  are  given  in 
Table  3.  As  may  be  noted,  the  use  of  the  moly  tab  results  in  a 
significant  reduction  in  chip  Q2fi  temperature.  The  moly  tab 
however  increases  the  temperature  of  chip  Q5B  because  of  the 
increased  thermal  conductance  between  Q2B  and  Q53*  Table  4  contains 
the  component  to  case  thermal  resistance.  The  significant  item 
shown  by  this  table  is  the  variation  in  the  thermal  resistance  of 
chip  Q5fi.  The  thermal  resistance  of  this  chip  is  directly 
influenced  by  the  power  dissipated  by  chip  Q2B*  Conduction  of  the 
thermal  energy  from  Q2:  to  Q5b  determines  the  resultant  temperature 
of  Q5b  and  thus  the  thermal  resistance.  As  such,  it  may  be 
surmized  that  the  use  of  thermal  resistance  is  not  always  an 
indication  of  a  component's  thermal  performance.  That  is, 
low  powered  chips  can  appear  to  have  a  high  thermal  resistance 
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FIGURE  9:  THERMAL  RESISTANCE  OF  SILICON  CHIPS  MOUNTED  ON  MOLT  TABS 


PARAMETERS  FOR  DESIGN  STUDY  OF  MULTIPLE  COMPONENT  MOUNTINGS 


MOLYTAB  DIMENSIONS  .065  X.  065  X  .005  INCH 
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CASE  1 

CASE? 

CASE  3 

CASE  4 

CHIP  <3® 

87.1 

87.7  •;> 

87.4  -■ 

40.9 

CHIPQjb 

4700 

mnjiu 

305 

9200 

RESISTOR 

513 

53.3 

53.7 

-  • 

— 

TABLE  4:  DESIGN  STUDY  COMPONENT  TO  CASE  THERMAL  RESISTANCE  °C/WATT 
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when  positioned  near  a  high  powered  chip,  while  still  operating  at 
acceptable  temperature  levels. 


3.4  Beryllia  Substrates 


A  current  state  of  the  art  technique  used  to  reduce  the  thermal 
resistance  for  high  powered  solid  state  elements  is  to  replace 
the  alumina  substrate  with  a  beryllia  substrate.  Beryllia  exhibits 
a  thermal  conductivity  nearly  equi valent  to  aluminum,  thus  the 
energy  dissipated  by  the  chip  is  diffused  through  out  the  substrate. 
The  results  for  a  beryllia  substrate  are  shown  in  Figure  11.  The 
chip  temperature  is  then  reduced  to  113. 5°C  in  comparison  to 
193°C,  for  an  alumina  substrate  with  the  same  boundary  conditions, 
Figure  12.  Beryllia  however  is  extremely  brittle  and  toxic  when 
being  worked.  Hence  it  would  be  desirable  to  find  another  means 
of  reducing  the  thermal  resistance. 


3.5  Alumina  Substrates  with  Integrated  Heat  Pipe  Cooling 


The  utility  of  using  heat  conducting  fluids  over  microcircuits 
was  found  to  be  ineffective  (Section  3.2).  It  was  then 
suggested  that  an  integrated  heat  pipe  cooling  concept  might  be 
able  to  reduce  the  chip  operating  temperature  and  decrease 
thermal  cycling  effects.  To  investigate  this  cooling  technique 
various  concepts  were  studied.  Figures  13  to  15  contain  illustra¬ 
tions  of  the  concepts,  thermo  physical  properties  and  results. 
Figure  13  depicts  a  single  chip  covered  by  a  5  mil  wick  saturated 
with  a  fluorocarbon  fluid.  The  chip  acts  as  the  evaporator  and 
the  surrounding  area  and  cover  form  the  condenser.  This  scheme 
reduces  the  chip  temperature  from  193°C  to  177. 5#C.  Thus,  this 
method  is  not  as  effective  as  the  moly  tab  mountings.  If  the  wick 
thermal  conductance  could  be  increased  by  an  order  of  magnitude 
then  the  chip  would  be  cooled  to  132.5*C.A  combined  moly  tab  heat 
pipe  concept.  Figure  15,  reduced  the  chip  temperature  to  152.7°C. 
Comparisons  between  this  concept  and  moly  tab  mounts  indicates  that 
the  integrated  heat  pipe  technique  does  not  offer  any  significant 
advantage . 
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FIGURE  11:  MICROCIRCUIT  WITH  BERYLLIA  SUBSTRATE 
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3.6  Alumina  Substrate  with  Sputtered  Beryllia 

Another  technique  for  diffusing  the  large  localized  heat  flux 
dissipated  within  a  solid  state  element  is  shown  in  Figure  16. 

This  concept  is  similar  to  moly  tab  mounts  wherein  the  beryllia  acts 
as  a  thermal  fin  to  reduce  the  heat  flux  incident  to  the  alumina 
substrate.  The  results  for  two  variations  of  sputtered  beryllia 
are  given  in  Figures  17  and  18.  The  chip  temperature  has  been 
reduced  significantly  and  this  concept  compares  favorably  with 
moly  tab  mounts  and  beryllia  substrates  techniques. 

3.7  Thermoelectric  Cooling 

It  has  been  proposed  that  for  Large  Scale  Integration  (LSI)  elect¬ 
ronic  modules  that  thermoelectric  coolers  be  applied  to  reduce 
solid  state  element  temperatures.  Figure  19.  Two  concepts  were 
analyzed  to  examine  this  cooling  technique,  a  single  chip  mounted 
on  an  alumina  substrate  and  one  using  a  beryllia  substrate.  The 
thermoelectric  geometric  configuration  and  materials  are  shown 
in  Figure  20.  Results  for  the  computation  are  given  in  Figures  21 
and  22.  The  chip  temperature  was  reduced  to  102.7  and  74.5®C 
respectively.  This  technique  offers  significant  advantages  in 
temperature  control,  however  the  COP  for  the  coolers  is  approximately 
(.371.  That  is,  for  each  watt  of  cooling  (2.691  watts  of  electrical 
energy  are  required  to  operate  the  coolers. 
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FIGURE  16:  BERYLLIA  SPUTTERED/THERMOELECTRIC  COOLED  ELEMENT 
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SCALE:  1  CM  ■  4  MIL 


FIGURE  18:  ISOTHERM  PLOT  FOR  MICROCIRCUIT  WITH  BERYLLIA 
SPUTTERED  OVER  TOP  SURFACE  AND  WITH  L  ■  30  Mi 
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FIGURE  20:  THERMOELECTRIC  COOLED  MICROCIRCUIT 
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FIGURE  22:  MICROCIRCUIT  WITH  THERMOELECTRIC  COOLER  -  BERYLLIA  SUBSTRATE 
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4 . 0  CONCLUDING  REMARKS 

The  numerical  results  described  in  the  previous  section  illustrates 
the  utility  of  thermal  models  in  predicting  component  temperatures 
and  in  aiding  the  engineer  during  initial  design  phases.  Various 
trade  studies  can  be  completed  rapidly  allowing  the  placement  of 
high  powered  components  in  regions  where  they  will  have  minimal 
effect  on  surrounding  devices.  Additionally,  different  mounting 
and  thermal  control  techniques  can  be  compared  in  an  effort  to 
chosse  those  offering  the  best  temperature  control. 

It  was  found  that  epoxy  "squeeze  out"  can  be  beneficial  in  reducing 
the  component  temperatures.  However  epoxv  squeeze  out  is  difficult 
to  control  and  thus  when  relied  upon  for  temperature  control  can 
cause  circuits  to  vary  widely  in  maximum  component  temperatures. 

The  use  of  molv  tabs  and  sputtered  bervllia  resulted  in  significant 
reductions  in  component  temperature  levels.  Additionallv,  it 
was  shown  that  the  integrated  heat  pipe  cooling  concept  does  not 
provide  effective  temperature  control. 

Thermoelectric  coolers  appear  to  offer  a  promising  technique  for 
controlling  LSI  module  temperatures.  This  method  however,  requires 
the  expenditures  of  significant  electrical  enerav  for  the  coolers 
to  operate. 
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APPENDIX 


Fundamentals  of  Thermoelectric  Coolj.n<? 


These  effects  are  associated  with  the  thermoelectric  effect: 
Seebeck 


Thomsen 


Peltier 

2 

Sometimes  Joulean  heating  Ci  R)  is  included  as  a  thermoelectric 
effect,  but  due  to  its  widespread  occurence  it  is  not  customary 
to  include  it  as  a  special  and  unique  phenomenon. 

The  Seebeck  effect  occurs  when 
2  conductors  are  connected  in 
a  circuit  with  their  junctions 

at  different  temperatures.  Higher  „cHfnchr 

temperatures  cause  more  electrons 
to  populate  states  above  the  Fermi 
energy  level.  Their  mobility  is  increased, 
and  a  diffusion  to  the  cold  end  occurs. 

With  more  electrons  at  the  cold  end  than 
at  the  hot  end,  an  electrical  potential 
is  set  up.  If  material  1  is  the  same  as  material  2,  the  same 
condition  is  set  up  in  each,  and  there  is  no  net  effect.  If  the 
materials  are  dissimilar,  the  excitation  to  higher  energy  levels 
is  different  in  each,  and  a  net  voltage  results.  This  effect 
gives  rise  to  the  "thermocouple  effect".  The  Seebeck  effect 
is  described  by  the  Seebeck  coefficient,  ol,  such  that 


oL  »  d dC,  _  r/e, 

'*  *r  dr  r*r  _ 

E12  is  called  the  Seebeck  potential  or  voltage.  1 1  00 

The  Thomson  effect  occurs  when  a  conductor 
is  subjected  to  a  thermal  gradient.  In  such 
a  conductor  electrons  would  migrate  to  the 
cold  end.  If  a  voltage  is  imposed,  the 
electrons  are  forced  to  move  in  the  direction 
of  the  potential.  In  the  case  where  electrons 

move  from  hot  to  cold,  electrons  give  up  enerqv  as  1 _ | 

heat.  If  the  electrons  are  forced  to  move  in  the 
opposite  direction,  i.e.,  from  cold  to  hot,  they  add  to  their 
thermal  energy  state  by  absorption  of  heat.  The  Thomson  effect 
is  described  by  the 
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Thomson  coefficient*  C  »  by 

O  Crtt—m*,*  )  *  Tl 

where  Q  is  the  rate  of  thermal  energy  absorption  or  emission,  I  is 
the  current,  and  dT/dX  is  the  temperature  gradient.  In  this  case 
Q  is  positive  when  added  for  4-1  in  the  +X  direction. 

The  Peltier  effect  is  the  reverse  of  the  Seebeck  effect 
in  the  sense  that  the  transformations  of 
energy  are  basically  caused  by  the  same 
phenomenon.  Electrons  convert  kinetic 
to  potential  energy  and  vice  versa,  and 
their  stored  potential  energy  depends 
on  these  environments,  i.e.,  the 
material  containing  them.  Therefore 
as  electrons  move  from  material 
1  to  2  they  take  up  or  release  energy 
according  to  the  capacity  of  the  com¬ 
bined  state  to  store  potential  energy. 

The  Peltier  effect  is  described  by  the 
Peltier  coefficient,  {or 

c'a^r  or  Q  *  X 


i 


i 

£ 


T 

Corrt*rf,  I 


Classical  thermodynamics  has  established  relations  among  the 
3  coefficients 

T 

1  O'  .  7-  &&•*) 

*  6/1  dT 

The  analysis  of  thermoelectric  cooling  systems  must  include  the 
irreversible  effects  of  heat  conduction  and  Joulean  heatinq. 
These  2  cam  be  described  by  the  relations 

Cone/ucArcn  Qc  =  (J  C  T;  -  71 ) 


where  U  is  an  overall  conduction  coefficient,  T.  refers  to  the 
hot  temperature,  and  Tq  refers  to  the  cold  or  source  temperature. 


Atm* ,  Oj  -  I  */Z. 

where  R  is  the  electrical  resistance. 

Consider  a  thermoelectric  system 
consisting  of  a  material  which  has 
a  positive  Seebeck  coefficient  (p)  and  a 
material  which  has  a  negative  coefficient 
Cn) .All  5  effects  occur  in  this  system: 
Peltier  cooling  Thompson  nesting 

heat  conduction  S’eebeck  voltage 

Joulean  heating 
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The  relation  f *  T  ^’'/e/r  shows  that  if  the  Seebeck  coefficient 
is  constant,  the  Thomson  coefficient  is  zero.  It  has  been 
shown  that  taking  an  average  value  of  cL  does  account  for  moderate 
Thomson  heating  effects. 

Also,  it  is  customary  to  assume  that  1/2  the  Joulean  heating 
goes  to  each  plate.  Therefore 

q.^u^t.1  -  u(r-r.)  -  ir'K. 

*  C,  *  cir.T;T-  ufr-T.)-  i  r‘it 


where  the  equivalence 
or  rearranging 

7~—  T  - 

'4  '  O 


r 

72 T  - 


has  been  used. 
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The  power  required  to  obtain  this  refrigerating  effect  is  due 
to  the  Seebeck  voltage  and  Joulean  heating 

W  -  (t-t.)]!  +  r^e 

The  coefficient  of  performance,  C.O.P.  is  defined  as  Q9 /W 


C.o.p 


olpo  7i  JO  -  0(77-71)  -  j 

oip*  (T'-T')r  * 


This  is  compared  to  a  Carnot  cycle  where  the  irreversible  effects 
are  absent  i.e.  if-  a  1  /?-  o 
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Optimizing  the  system  for  Q  or  T4-T  with  respect  to  current 
results  in  o  i  o 


- 


Also  if  we  set  Q  ■  0  for  the  maximum  CT.-T  )  and  use  the 
optimum  current 


r--r**v(cL^T.I-{T'Z.)  3  £7  ( n 
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The  "figure  of  merit"  is  defined  as 
such  that 


V < 


Z  is  a  function  of  material  properties  and  typical  values  are 
given  in  the  table  belov 

Material  2  1/°K 


Chromel-constantan 

0.0001 

Sb-Bi 

.00018 

Sb-C91%  Bi,  9%  Sbl 

.00023 

ZnSb  -  Constantan 

.0005 

PbTe  (p  type!  -  PbTe  Cntype) 

.0013 

Bi-Te3  (pI-Bi-  Te3  Oil 

.002 

p-n  Alloy  Thermocouple 

C.50%  Bi2Te3,  50%  Sb2Te3J 

.0025 

C75%  Bi  Te3,  25%  Bi2Se3 

+  minor  CuBr2  added! 

.0025 

For  the  above  values  of  Z,  0  *  k,  thermal  conductivity,  and  R  *  , 

electrical  resistivity. 

3  C  ^  ^  o 

To  optimize  C.O.P.  with  respect  to  current,  jf 

CoLT.-xtzX*  atx+z'x.)  -  (,tT.x  - ua r- iz te Ydur X* )  ^ 

<=>x  "  C  *>.  (r..r„)x*  r*«  ;v 

multiplying  out  and  collecting  terms 

xx(/z>xr»  +  f  at*.  )-  fi  u  **r)x  -it*  zru  =  o 


After  some  arithmetic ,  solving  the  quandrtic,  and  using  the 
definition  of  Z 


ci*T(t  *  Yi+  ZX*’) 

Z  ' 


7T*-7Z 

z 


The  +  sign  applies  _  _  dAT  (/  f-rT^TTL  Yj-  /77T7Z) 

and  multiply  O-yKJK)  *  ~  G2TZ.  — ft-ft+l'K) 


***  = 


dAT  / 

(YTT7?, r-,) 


Urlvtk.  7J.  * 
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Straightforward  substitution  of  t w  into  the  above  equation  for 
C.O.P.  yields  the  maximum  C.O.P.  ,.s.  e.o./t.  Ar  jr^ 


-~r  f/TTiJZ  -  f ) 

XL  \  _  K-7j  ( ' _ ’»/ 


+  / 


cc.o.p. Imax  and  Cr^~T0Ln|ax  are  plotted  on  the  next  page  for 
T.  »  300°k  and  various  values  of  S. 


Summary: 


AT 


x(y77Fn. r  -/) 


V-r 

R 


(co.fi.) 


nmfi 


Sr  -  f.  ) 

/TTTnT  *•  * 


max.  o.  o.  p. 


4m  S’ 


-=  jr  tT  %  (*o  /o*d ) 


’/cr  max. .  7“ 

^  max.  Q 


Note  that  an  estimate  of  z  can  be  made  by  running  at  "no  load" 
and  measuring  T  &  T  .  The  above  expressions  are  for  a  specific 
value  of  Z  *  /R.u  ,  but  R  and  U  involve  geometry  which  may  be 
optimized  for  performance.  Consider  U,  the  overall 

conduction  path 

u~  K  £•  +  kr  &  *  K  + 

uJitrt.  4  «  <ti tm.  t  L  •  /tfft*.,  j  *  fi  /L 

Rm  fn  /?„  +•  f*  /fy 

■SolcJi  Ru*  fa  -  fa  (ft /fy )  +■  &rfn  ((>/?*)*  fi* 

7o  c  RU  tail.  *-q  t  a  net  -Hi+i- 

f1*  /fp  *  ( f*  ^  / 1 Pp  ¥-A<*  acne*. *><rx.  .  Rsr  c5**  /<J^ 

Am  /jC*  * 
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Maximum  theoretical  temperature  difference  for 
thermoelectric  element. 


Temperature  difference.  'K 


The  maximum  coefficient  of  perform¬ 
ance  that  can  be  achieved  with  a  given  temperature 
difference  and  a  specified  figure  of  merit  for  the 
couple.  The  hot  shoe  temperature  is  assumed  to  be 


Maximum  theoretical  coefficient  of  performance  for  a 
thermoelectric  element. 


If  an  electrical  contact  resistance  of  r  is  assumed  at  each  junction 
Cr  at  each  end  of  each  element 1 


Also  for  RU  to  be  a  minimum, 

*°m  l (&*”)*  +-  (fi-krf  ] 

If  we  had  2  materials  with  the  same  aL  ,  k,  and  p  ,  we  could 
write  J 


Some  of  the  above  properties  are  governed  by  the  Wiedermann- 
Franz  Law  kf  *  CLorentz  number)  T.  For  metals,  the  Lorentz 

—8  2 

number  has  a  value  of  2.5  x  10  watt-ohm/°K  .  Also  for  metals, 
the  Seebeck  coefficient,#,  has  values  on  the  order  of  10~5  volt/ 
°C.  Therefore  ^  .-cL*/pk  ir  (/o~*) /fas’* /o~*)  r  x.  **/©'  >/t. 

For  temperatures  on  the  order  of  300°K,  2  10-5.  It  is  apparent 
that  thermoelectric  cooling  with  metal  conductors  is  of  minimal 
value . 

Values  of  properties  are  shown  below  for  some  metals  and  semi¬ 
conductors  where  ora  p_/ ,  the  electrical  conductance . 


TABLE  .2 

Thermoelectric  properites  of  metal 


Electrical  Seebeck 

Thermal 

« 

Figure  of 

Conductivity  Coefficient 

Conductivity 

cL  T 

Merit 

a. 

fiu+H-  /Cm  'C  ) 

copper  5.2  X  lof 

+2.5 

3.96 

3.7  x 

10"| 

93  x  10 

Nickel  1.5  x  IO3 

-18 

0.87 

4 . 9  x 

10-5 

5  6  x  10 

Bismuth  8.6  x  10 

-75 

0.08 

4 . 8  x 

10  5 

6  x  10 

TABLE  3 

Thermoelectric  properties  of  semiconductors 

<r  * 

,  j*.  fc' 

(L %-C+m.)  'S’ 

<  u 

kM*  + 

d  \ 

r 

,  ¥L  V 
a*  ts- 

*. 

*  "1C 

c* *»  *e 

Germanium 

1000  200 

0  63  0  006 

4  x 

10-5 

f  3  x  10 

Silicon 

500  200 

1  13  0  003 

2  x 

10_5 

1  8  x  10 

Indium  antimonide 

2000  200 

0  16  0  01 

8  x 

iO-4 

4  7  x  10 

Indium  arsenide 

3000  200 

03  0  015 

1  2  x 

10-5 

3  8  x  10 

Bismuth  telluride 

1000  220 

0  016  0  004 

4  7  x 

10  5 

2  3  x  10 

2T»  r#8LE.  3  ,  4L  ^  /•#*% 

Oenetuc¥*s'/y  ¥ 

k  •  4**4 

* 

4 
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Insulators  Semiconductors 


Metals 


A  schematic  representation  of  how  the  F*ebeck 
coefficient,  resistivity,  and  thermal  conductivity  depend  on  the 
concentration  of  extrinsic  charge  carriers. 


The  variation  of  a£,p,  k,  and  Z  as  a  function  of  extrinsic  charge  carrier 
is  sketched  above.  kg  is  the  electronic  portion  of  the  thermal 

conductivity,  and  k,  is  that  part  associated  with  the  lattice  or 
phonon  transfer. 

The  table  below  is  from  Spring,  but  note  that ^  can  not  be  separated 
out  as  an  independent  quantity. 


Thermoelectric  materials 


-Spring: Direct  Generation  of 


1 

2 

3 

4 

5 

0 

7 

8 

Electricity 

9  10 

(eV) 

Melting 

point 

CO 

Type  of 
semi* 
conductor 

a?o* 

K* 

(W/cra  deg  C) 

Ki 

(W/cm  deg  C) 

Tempera¬ 

ture 

(CK> 

-Maximum 

operating 

temperature 

(°K> 

BIjTej 

015 

575 

n  or  p 

4  X  10-* 

0  004 

0016 

2  X  10-* 

300 

450 

BiSb.Te7l| 

— 

— 

P 

4-6  x  10- 5 

0  004 

0010 

3-3  X  10- * 

300 

450 

BijTejSe 

0-3 

— 

n 

3-6  x  10*  * 

0003 

0013 

2-3  x  10*  * 

300 

000 

PbTe 

0-3 

m 

nor  p 

2-8  X  10* s 

0003 

008 

1-2  X  10-* 

300 

POO 

GeTe<+  Bi) 

— 

725 

P 

3-2  x  10- » 

— 

002 

1-6  x  10-* 

soo 

900 

ZnSb 

0-6 

546 

P 

2  X  10" * 

— 

0-17 

1-2  X  10-* 

500 

600 

AgSbTea 

oa 

576 

P 

1-4  X  10*» 

0002 

0006 

1-8  X  10'» 

700 

900 

InAaf-t-  P) 

0-45 

040 

n 

4  X  10-» 

— 

0-07 

6  x  10-* 

900 

1100 

CeS{+B») 

— 

— 

n 

8  X  10-* 

— 

001 

8  x  10- ♦ 

1200 

1300 

Cu,TejS 

— 

030 

— 

18  X  10*« 

— 

0012 

1-5  X  10- * 

1100 

Ge-Si 

0-8 

— 

n 

P 

/  3  X  10*‘ 
\  2  x  10-» 

— 

0033  ^ 

'  9  X  10‘« 

6  X  10-* 

900 

900 

1200 

The  figures  given  in  columns  5,  8,  7  and  I  refer  to  the  temperature  quoted  in  column  9. 
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2(° K)-'  X  |03  p(ohm-cm)  X  I03  <»(yV/°K) 


P 


* 


.M  I  4  ■  I  |l|  .  ■  I 


/ 


Temperature  (°C) 


The  Seebeck  coefficient,  electrical  resistivity, 
and  figure  of  merit  of  a  p-type  alloy:  AgSbTe;. 


The  Seebeck  coefficient,  electrical  resistivity,  and 
figure  of  merit  of  an  /(-type  alloy :  75  %  PbTe,  25  %  SnTe. 


I 
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performance  as  previously  defined: 


2 


(T~) 


In  the  optimization  we  use  the  Wiedermann-Franz  Law 
where  L  =  Lorenz  constant  C2.5  x  10” 9  /  •k'-  A*  ' 

Therefore  4 

f  — — I  =.  cons+a»+  (\m/.d  Jcr  fr^-T,)  i  ( o 

I  ft*  Yl*  J 


and  to  optimize  ju. 


It  is  convenient  to  use  Y  *  d  8=1 
with  respect  to  d.  If  this  is  done, 

(i\-  [&),(%),]*-  (M 

Also  it  can  be  determined  that  for  max.  c.o.p. 


fax 


V* 


(^)i  Yv-y-.) 

Ki  (S'  -•) 
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With  respect  to  temperature  distribution, /a  is  a  minimum 

when  *  - -  77., 


The  above  optimization  procedures  make  the  c.o.p.*s  for  all 
the  stages  equal.  For  these  conditions 


The  proceeding  theory  was  coded  and  used  as  a  subroutine 
to  MZCTA.  A  listing  is  given  below. 


1 


LEVEL 

21.6  I 

NAV  72 

»  C./363  FC*T*AN  N 

compile* 

OPTIONS  -  NAME*  MA  IN*  OP  T*02 . L INECNT»54  .  S IZE* jOOOK , 

S OUPCE. EBCDIC. ':a  rr, NOO^m.LCAO.MA'’, NOEDIT.  rO.NOXREF 

ISN 

0332 

SUBROUTINE  COOLER  ICN.T) 

ISN 

3003 

DIMENSION  fNIU.TIl) 

ISM 

000 + 

C 

REAL . I3PT 

A  IS  THE  SEE8ECK  C06F.  tN  VOLTS  PER  0 EG  KELVIN 

ISN 

OOOS 

c 

A  *  .30038 

8  IS  THE  PESISTIVITV  OF  BI2-TE3  I N  'JHM-MILL 

ISM 

OCOfc 

9*. 3937 

ISN 

0007 

APN  *  CNl  9)  *CN(  24 1 

ISM 

3338 

R  -  2.*8*ICNI75I *CN 1 761 )/ APN 

ISN 

0009 

IOPT  *  8.0 

ISN 

0010 

Of 03  »  A*ICPT*ltl 12691*273. 1-.S*I0PT**2*R 

ISN 

0011 

•0004  ■  A*10PT*I  71  12721*273. >-.5*ICPr**2*R 

ISN 

0012 

0C39  *  A*IOPT*|T| 13291*273. I-.5*ICPT**2*P 

ISN 

3313 

Of 13  *  4*I0PT*IT< 13321*273. 

ISN 

0014 

• 

OH) 3  •  A*I0PT*|TI1469)*273. »*.5*ICPT**2*R 

ISM 

001S 

QM04  •  A* IOPT*  1  Tl  14 721*273.  1  *.5* I QPT*1*2*P 

ISN 

3016 

3M09  •  A*IOPT*ltl 15291*273. )*.5*lOPT**2** 

ISN 

0017 

OHIO  »  A*ICPT*ITI 1532 1*273. I*.5*IUPT<>*2*» 

ISN 

0018 

CNl 4731 )»QC33 

ISN 

0019 

CNl 4792 1*0004 

ISN 

0320 

CNl  47031*0009 

ISN 

0021 

CM47341-OCIO 

ISN 

0022 

CNl 4705 1*0H03 

ISN 

0023. 

CNl 4736 l*OM04 

ISN 

0324 

CNI473TI*0M39 

ISN 

0025 

CN| 4702 )*0H10 

ISN 

0026 

CM| 4751  l»R 

_ ISM  3027 _ 

CNI4752 l-APN 

ISN 

0028 

CNl  475i  I*  IOPT  I 

ISM 

0029 

RETURN 

ISN 

0030 

END 
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